#pragma once

void print_matrix(double** phi, int n, char* buf);

void print_matrix2(double** phi, double** eta, int n, char* buf);

void print_matrix3(double** phi, double** eta, double **xx, int n, char* buf);

void print_vector(double *x, int n, char* buf);

void Cholesky(double **a, double **L, int n);

void equ_solver(double **a, double *x, double *b, int n);

void Cholesky_optimal(double **a, double **L, double *d, int n);

